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(54) Method for detecting rare events. 



(57) This invention relates to a method for the 
detection of rare events which are events occur- 
ring at a frequency of less than one in 10 4 . It is 
particularly useful in detecting events at fre- 
quencies of one in 10 5 and one in 10 6 . The 
method employs one or more first markers that 
are specific for the rare event particle each of 
which are labelled with a dye having an emis- 
sion wavelength distinguishable from the 
other(s). The method further employs one or 
more second markers that are specific for a 
majority of the remainder of the particles pre- 
sent in the cell sample but are negative for the 
rare event particle. The second markers all are 
labelled with the same dye. The second markers 
are collectively referred to as the exclusion 
color. By analyzing the particles for the pre- 
sence of the first marker(s) and the absence of 
the second markers, rare events can be detec- 
ted. 
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Field of the Invention 

This invention relates to a method for the detec- 
tion of rare events in a sample of particles, and more 
particularly relates to a method for the detection of 
rare cells in a sample of cells at frequencies of less 
than one in 10 6 cells. 

Background of the Invention 

Particle analysis generally comprises the analy- 
sis of cells, nuclei, chromosomes and other particles 
for the purpose of identifying the particles as mem- 
bers of different populations and/or sorting the par- 
ticles into different populations. This type of analysis 
includes automated analysis by means such as flow 
cytometry. In flow cytometry, the particle, such as a 
cell, may be labeled with one or more markers and 
then examined for the presence or absence of such 
markers. In the case of a cell, such as a leukocyte, tu- 
mor cell or microorganism, the marker can be directed 
to molecules on the cell surface, such as receptor pro- 
teins, or to molecules in the cytoplasm or nucleus 
such as proteins, enzymes or nucleic acids. Examin- 
ation of a particle's physical characteristics ( e.g. , 
size), as well as for the presence or absence of mark- 
ers) can provide additional information which can be 
useful in identifying the population to which a particle 
belongs. 

Flow cytometry comprises a well known method- 
ology using multi-parameter data for identifying and 
distinguishing between different particle types in a 
sample. For example, a sample of cellular particles 
may be drawn from a variety of biological fluids such 
as blood, lymph, bone marrow, ascites, pericardial 
fluid, pleural fluid, cerebrospinal fluid or urine, or may 
be derived from suspensions or aspirates of cells 
from organs and solid tissues such as brain, Gl tract, 
prostate, lung, breast, kidney, lymph node or liver. In 
a flow cytometer, cells are passed substantially one 
at a time through one or more sensing regions where 
in each region each cell is illuminated by an energy 
source. The energy source generally comprises an il- 
lumination means that emits light of a single wave- 
length such as that provided by a laser ( e.g. , He/Ne 
or argon ion) or a mercury arc lamp with appropriate 
bandpass filters. Light at488nm is a presently prefer- 
red wavelength of excitation when a single light 
source is used and 632nm is preferred as a second 
light source when two light sources are used. 

In series with one or more sensing regions, mul- 
tiple light collection means, such as photomultiplier 
tubes, are used to record tight that passes through 
each cell (generally referred to as forward light scat- 
ter), light that is reflected orthogonal to the direction 
of flow of the cells through the sensing region (gen- 
erally referred to as orthogonal or side light scatter) 
and fluorescent light that may be emitted from the 



cell, if it is labeled with fluorescent marker(s), as it 
passes through the sensing region and is illuminated 
by the energy source. Each of forward light scatter (or 
FSC), orthogonal light scatter (SSC), and fluores- 

5 cence emissions (FL1 , FL2, etc.) comprise a separate 
parameterfor each cell (or each "event"). Thus, for ex- 
ample, up to four parameters can be collected (and re- 
corded) from a cell labeled with two different fluores- 
cence markers ( i.e. , two scatter parameters and two 

10 fluorescent parameters per cell). 

Flow cytometers further comprise data acquisi- 
tion, analysis and storage means such as a computer 
equipped with appropriate software. Signals from 
each PMT are converted into multiple data channels 

15 and stored by the computer. The purpose of the ana- 
lysis system is to classify and count cells wherein 
each cell presents itself as a set of digitized parame- 
ter values. Typically, by current analysis methods, 
data collected in real time (or recorded for later ana- 

20 lysis) is plotted in 2-D space for ease of visualization. 
Such plots are referred to as dot plots and a typical 
example of a dot plot drawn from light scatter data re- 
corded for leukocytes is shown in FIG. 1 of U.S. Pat. 
No. 4,987,086. By plotting orthogonal light scatter 

25 versus forward light scatter, one can distinguish be- 
tween granulocytes, monocytes and lymphocytes in a 
population of leukocytes isolated from whole blood. 
By electronically (or manually) gating on only lympho- 
cytes using light scatter, for example, and by the use 

30 of the appropriate monoclonal antibodies labeled with 
fluorochromes of different emission wavelength, one 
can further distinguish between cell types within the 
lymphocyte population ( e.g. , between T helper cells 
and T cytotoxic cells). U.S. Pat. Nos. 4,727,020, 

35 4,704,891 , 4,599,307 and 4,987,086 describe the ar- 
rangement of the various components that comprise 
a flow cytometer, the general principles of use and 
one approach to gating on cells in order to discrimin- 
ate between populations. 

40 While the above-described methods are useful 

and accurate for the analysis of particles that occur 
in a sample at a frequency of greater than 1%, accu- 
racy is much less reliable in the analysis of particles 
that occur with lesser frequency. These events are 

45 termed "rare events." 

Accuracy in the counting and analysis of rare 
events is important. For example, in many cancers 
elimination of "cancer cells" is the desired goal of a 
therapy. It is important to know if all of the cancer cells 

so have been eliminated. If a number of cancer cells re- 
main but are undetected because existing instrumen- 
tation cannot reliably detect the cells that remain be- 
cause of their low frequency, a patient could be given 
an erroneous "clean bill of health." Similarly, if these 

55 minimum residual cells have not been eliminated, it is 
important for monitoring purposes to detect increases 
in the frequency of such cells at the earliest possible 
time in order to determine whether changes in treat- 
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ment are required. Thus, accurately quantitating the 
number of cancer cells that remain after treatment is 
important in the future treatment and monitoring of 
cancer. 

A concern in the accurate counting of cells by 
flow cytometry is assuring that an event which the in- 
strument senses as being positive for a "fluorescent" 
marker actually is so as a result of being labelled with 
a fluorescent marker. It has been reported that fresh- 
ly isolated unlabelled mouse bone marrow contains 
autofluorescent particles that a cytometer senses as 
being "fluorescent." The incidence of these events 
was reported as one per 10 4 cells. As a result, it has 
been considered that the lowest achievable level of 
detection using commercially available cytometers is 
about one cell in 10 6 . 

In addition to unlabelled cells being recorded at 
some frequency as "fluorescent" events, the use of 
immunofluorescence markers can cause additional 
errors in data analysis. Monoclonal antibodies label- 
led with fluorescent dyes have a high affinity for their 
antigen. Unfortunately, because of the nature and 
structure of the cells on which such antigens appear, 
immunofluorescence markers will bind to structures 
other than the intended antigen. Thus, irrelevant cells 
can be recorded as positive fluorescent events due to 
non-specific labelling. For cells whose frequency is 
below 1% of the cells in a sample, this is an additional 
problem to be addressed. 

A third problem with conventional cytometry is 
the speed at which cells are analyzed. Current instru- 
mentation operates at an acquisition rate of 5,000 to 
10,000 events per second. Current instrumentation 
also has a limited random access memory associated 
with the data analysis hardware to acquire events for 
real time analysis. Thus, to acquire about 10 8 events, 
it takes about 3 to 6 hours at a minimum. To maintain 
cells over such a long period of time with the same flu- 
orescent characteristics, it is necessary to fix the 
cells. Fixation will have some effect on the physical 
characteristics of the cell as opposed to fresh living 
cells, and thus, presents a further opportunity for the 
introduction of error into the system. 

Accordingly, an improved method for the analysis 
of cells is required if rare event analysis is to be a 
more useful and accurate tool in cytometry. 

Summary of the Invention 

This invention comprises a high-speed method 
for the detection of particles in a sample at frequen- 
cies of less than one in 10 4 . 

The particles can be cells in biological fluids, in- 
cluding red and white blood cells, dis-aggregated 
cells from organs and solid tissues and microorgan- 
isms. Samples of such cells can be taken from blood, 
urine, bone marrow, ascites, pericardial fluid, pleural 
fluid, cerebrospinal fluid, lymph and other biological 



fluids, as well as from organs and solid tissues such 
as brain, kidney, prostate, liver, spl en, lung, breast, 
lymph node and Gl tract. 

The particles to be analyzed are combined with 
5 at least one first marker that is specific for ( i.e. , has 
a high binding affinity for) the particles to be detected 
and at least one second marker that is specific for the 
majority of the remainder of the particles in the sam- 
ple but is negative for ( i.e. , has low or no binding af- 
10 f inity for) the particles to be detected. It should be ap- 
preciated that in few circumstances will there be a sin- 
gle marker which will be specific for the majority of re- 
maining particles in a sample. Therefore, most uses 
will require more than one second marker. This sec- 
ts ond marker also can be referred to as the exclusion 
marker. 

Where the particles comprise cells, it is preferred 
that the markers be immunofluorescence markers. 
Where more than one immunofluorescence marker is 

20 used, it is preferred that each of the first markers be 
labelled with a different -fluorescent dye. Each such 
dye should have an emission wavelength that is dis- 
tinguishable from the other(s). It further is preferred 
that the second marker be labelled with a fluorescent 

25 dye that has an emission wave length distinguishable 
from that of the dyes used in connection with the first 
marker(s). It should be noted, however, that while 
each of the first markers generally is labelled with a 
different dye (when multiple first markers are used), 

30 when more than one second marker is used they all 
preferably are labelled with the same dye. Finally, it 
is to be appreciated that where multiple markers are 
used they can comprise combinations of immuno- 
fluorescence markers and other cytoplasmicior nu- 

35 clear markers ( e.g., nucleic acid dyes). - 

In the preferred embodiment, the number of first 
markers is two, and the number of second markers 
is greater than one. 

This invention has particularly utility in the de- 

40 tection of cancer cells in bone marrow or peripheral 
blood preparations of stem cells. Stem cells have 
been described as being CD34 + . Subsets of such 
cells include cells that are CD34 + /CD38~ cells that 
are CD347CD38"/HLA-DR + and cells that are 

45 CD34 + /CD38"/HLA-DR". In the treatment of patients 
requiring infusion of stem cells, it is important to elim- 
inate any remaining cancer cells from the prepara- 
tion. For those cancers that are CD34~, antibodies 
against CD34 would comprise the exclusion color 

so while antibodies specific for the cancer cells would 
comprise the first marker(s). 

Brief Description of the Drawings 

55 FIG. 1 is a plot of natural log probability of each 

sequence of 1000 events for a sample of 140 x 10 6 
cells containing rare events at a frequency of about 
10- 6 . Each of the probabilities is calculated using the 
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relative frequencies of parameter values in accor- 
dance with the algorithm set forth below. 

FIG. 2 is a schematic diagram of the gates ap- 
plied to n-dimensional data (and shown in 2-D form) 
in order to isolate rare events. 

FIG. 3 comprises a series of dot plots of (A) - (G) 
SSC versus FSC, (H) FL2 versus FL1 and (I) FL2 ver- 
sus FL3 for various amounts of REH cells seeded into 
10 s PBMC. 

Detailed Description of the Invention 

Peripheral blood mononuclear cells ("PBMC") 
were isolated from buffy coats by density gradient 
separation. After two washes with 0.45p, filtered Ca^ 
and Mg 2+ free phosphate buffered saline ("PBS"), the 
cells were resuspended in 0.45n filtered PBS in 20% 
fetal calf serum ("FCS"). Cell suspensions of more 
than 1.2 x 10 9 cells were rotated for one hour. Viability 
was determined by the propidium-iodide exclusion 
test on a FACScan brand flow cytometer (Becton 
Dickinson Immunocytometry Systems, "BDIS"). Vi- 
ability was always greater than 98%. 

Cells not stained were counted on a hemocytom- 
eter. 2.5 x 10 8 PBMC were aliquoted in up to 4 sepa- 
rate 15ml tubes. Each tube was supplemented 0.45u. 
filtered FCS to a total volume of 14ml. 

The pre-B cell line REH was cultured in RPMI in 
10% FCS in an atmosphere of 7% C0 2 at 37°C. 
Twelve hours prior to use, the cultures were split and 
fresh medium was added. The cells were harvested 
and washed with 0.45|i filtered PBS prior to sort. Vi- 
ability was determined as above. Viability was al- 
ways greater than 94%. 

Where required, pre-B cells were labelled with 
immunofluorescence markers as follows. Pre-B cells 
were sorted as described below into 2.5 x 10 8 PBMC. 
The mixture then was rotated for one hour to ensure 
dispersion. The cells then were spun down and su- 
pernatant removed to leave about 50u.l of solution. 
50|_lI of normal mouse serum (Dako) was added. The 
cells were resuspended and rotated for one hour. 

190(il of labelling cocktail was added. The cock- 
tail contained 25% normal mouse serum and the fol- 
lowing antibodies at the following final sample con- 
centrations. To separate pre-B cells from mature B 
cells, anti-CD19 fluorescein isothiocyanate (Leu 12 
FITC, BDIS) 4.6ng/ml and anti-CD20 FITC (Leu 16 
FITC, BDIS) 2.4u.g/ml were added as a first marker. 
To exclude T cells, granulocytes, monocytes and pla- 
telets, anti-CD7 peridinin chlorophyll protein (Leu 9 
PerCp, BDIS) 22.9|ag/ml, anti-CD3 PerCp (Leu 4 
PerCp, BDIS) 5.4u.g/ml, anti-CD4 PerCp (Leu 3a 
PerCp, BDIS) 1.6u.g/ml, anti-CD8 PerCp (Leu 2a 
PerCp, BDIS) 2.0u.g/ml, anti-CD15 PerCp (Leu M1 
PerCp, BDIS) 4.7ug/ml, anti-CD15 PerCp (L-16 
PerCp) 6.4 ug/ml, anti-CD14 PerCp (Leu M3 PerCp, 
BDIS) 1.3n g/ml, anti-CD41a PerCp 12.3u.g/ml and 



anti-CD42b PerCp 29.5m- g/ml were added as second 
markers. For additional identification of pre-B cells, 
anti-CD10 phycoerythrin (anti-CALLA PE, BDIS) 6.2 
jag /ml was added as an additional first marker. As an 

5 isotype control, anti-KLH PE 6.2u.g/ml was added. 

The cells then were resuspended and rotated for 
one hour to ensure complete labelling. The cells were 
washed twice in PBS and fixed in a total volume of 
200ml in 0.5% paraformaldehyde for 12 hours. After 

10 harvesting, the cells were washed once in PBS and 
finally resuspended in a total volume of 45ml 
PBS/0.25% FCS. The yield of recovered PBMC was 
58.2 ± 12.5%. 

All cells were analyzed on either a FACScan 

15 brand flow cytometer or a FACstarP'" 3 brand cell sor- 
ter (BDIS). The instruments were equipped with 
488nm argon lasers. The cell sorter was further 
equipped with a 70n nozzle with a drop drive frequen- 
cy of 27 kHz with 3 drop envelopes. 

20 Data acquisition and analysis was performed on 

a Hewlett Packard 340 computer with 16Mb of RAM 
using Lysis software (BDIS). Pulse width measure- 
ments were used to discriminate single cells from cell 
aggregates. 

25 Prior to any experiment being run, each instru- 

ment was cleaned to assure removal of particles and 
debris from prior experiments. This process is descri- 
bed in U.S. Patent No. 5,087,295. Briefly, a mixture 
of 0.07M NaOCI and 0.05M NaOH in particle free 

30 deionized water were added to the sheath tank. The 
instrument then was run for 20 minutes by-passing 
the built-in 0.2jj, sheath fluid filter. This tank then was 
emptied and particle free deionized water was added 
and the instrument then run for 30 minutes. The tank 

35 was again emptied and sheath fluid then was added 
to the tank. The fluidics system was allowed to equi- 
librate for 1 5 minutes thereafter. 

Electronic stability of the FACScan instrument 
was checked by a continuous run for 102 hours. The 

40 instrument has a test mode which sends sequenced 
light pulses (approximately 2000/second) to all PMTs. 
Alive gate was set in FL1 and FL2 in accordance with 
manufacturer's instructions. This gate excluded all 
events at or below the light pulse intensities and indi- 

45 cates spurious events generated by the electronics of 
the instrument. No spurious events were generated 
in the test run. 

For rare event acquisition, approximately 10 8 
PBMC were run on an instrument at an acquisition 

so rate of 5,000 to 6,000 events per second. Event ac- 
quisition was triggered on a threshold set in FL2. The 
threshold was set just above the upper limit of auto- 
fluorescence of the bulk of the negative cells. A 
strong FL2 signal on possible rare events was prefer- 

55 red to increase the probability of detection of rare 
events. It should be appreciated that fluorescence 
thresholds may be set for other fluorescence parame- 
ters as well in order to further refine the selection 
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method. In addition, a live gate was set in forward and 
side scatter which rejects most of the very small par- 
ticles such as platelets. These procedures reduce the 
number of events taken to RAM and to the system's 
hard disc to about 1%. On average, 10 6 events were 
recorded for about 10 8 cells analyzed. 

Prior to the analysis of the acquired data file, an 
algorithm was applied to remove bursts of events. 
FIG. 1 displays event number versus the natural log 
of the probability of each sequence of 1 000 events for 
a sample of 140 x 10 6 cells containing rare cells at a 
frequency of 10" 6 . The black heavy lines at the top of 
the figure represent rare events acquired as defined 
by the gates described in FIG. 2. 

Three major bursts are apparent. The bursts of 
events were not of a Poisson distribution, and thus, 
represented clumps of particles or debris the inclu- 
sion of which would lead to spurious results. In two of 
the bursts, "rare" events were detected. The algo- 
rithm described below removes all events inside of a 
burst, and thus, excludes those "rare" events from 
being taken into account. The conjecture is that these 
events are not real rare events, or might not be, and 
therefore should be disregarded to assure integrity of 
the system. It should be readily apparent that this al- 
gorithm has special utility for rare event detection but 
will have general utility to the acquisition and analy- 
sis of any data set where fidelity is desired. 

The algorithm was applied as follows. Assume a 
data set with L events consisting of N parameters 
digitized in B bits. Then each event Eg (i=1 . . L) will 
receive a probability Pj based on a number of W 
events in a window around Ej. Let H tjk be the number 
of events in the window around Ej in which the jth 
components have a value k (i=1...L, j=1.. .N, 
k=1...2 B ). Then H ljk is Poisson distributed with X lk - 
=WP jk , where P jk is the a priori probability that para- 
meter j has value k. This applies only if X jk is moder- 
ate. Pj k can be known or estimated. In the later case, 
it can be estimated during a test run or from the data 
itself as a ratio of the total number of events with the 
kth component equal to j and the total number of 
events L. 

Any realization of H ijk has a certain probability 
associated with it, namely: 



Pijk = e( - V 



Hijk 1 



Multiplying the P ijk over j and k gives the desired 
probability P iP if all P ijk were independent (which they 
are not because the number of events in a window is 
limited). However, if a fraction f of P ijk with f « W is 
taken, then all of the probabilities can be taken as in- 
dependent when the f is spread equally over all pos- 
sible j. 

Once all Pi have been found, a minimum filter of 
size W is applied. The minimum filter assigns to each 
event the minimum value of all the P, in the window 



around it. This has the effect of eliminating certain 
events that occur just before or after a burst. 

Finally, a quantile Q of all Pj is found. (Ten per- 
cent provides an appropriate measure for Q.) All 

5 events with a P| less than Q are removed from the 
data set resulting in a data set of length (1-Q)*L. It 
should be noted that the algorithm will still remove Q 
events; however, because the probabilities for these 
events vary only due to statistical noise, there will be 

10 no consequence for data analysis. 

Modifications of the method to calculate the 
probabilities ( e.g. , considering each event in n-di- 
mensional space) or a direct "is in agreement with" or 
"is not in agreement with" the distribution as deter- 

15 mined from the whole data file. Furthermore, one 
could use a different method to find the threshold be- 
low which events are removed ( e.g. , visual inspection 
of the calculated probabilities). 

Referring to FIG. 2, as general outline of how first 

20 and second markers are used to identify rare events, 
such as white blood cells, is shown. First, a fluores- 
cence threshold is applied above the autofluores- 
cence background of the cells to exclude,all events 
that do not have an immunofluorescence marker. 

25 Second, a region, R 1p is found in FSC and SSC that 
will exclude small particles, such as platelets. See 
FIG. 2A. It should be appreciated that if the exclusion 
of small particles is not necessary or may accom- 
plished through the use of the second markers, this 

30 step can be omitted. Finally, fluorescence regions, R 3 
and R 5 , are found for at least two combinations of the 
fluorescence data where rare events should occur. In 
this example, for pre-B cells, the plot is of the two first 
markers ( i.e. , FL2 versus FL1 ) which will contain only 

35 labelled B cells (see FIG. 2B), and there is a plot of 
one first marker versus the exclusion marker(s) ( i.e. , 
FL2 versus FL3) which will contain only labelled B 
cells (see FIG. 2C). A rare event, thus, is defined as 
a cell exceeding the fluorescence threshold and fall- 

40 ing within each of the 3 gates set. 

It should be apparent to one of ordinary skill that 
these displays are not necessary for data analysis, 
but are provided as a means to visualize the steps in- 
cluded in the method of the invention. These steps 

45 generally are run without the need for the user to see 
them. It also should be apparent that as the number 
of first markers increases, so too will the number of 
gates set for the various combinations of first mark- 
ers. 

so The following examples described in more detail 

the nature of this invention. 

Example I 

55 250, 2,500 and 25,000 pre-B cells were sorted 

into 2.5 x 10 8 PBMC. This represents a concentration 
of pre-B cells of 10 10- 5 and 10- 4 respectively. 
Sorting was done using a FACStar p,us cell sorter us- 
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ing the instrument's "count mode" and "pulse proces- 
sor." Count mode delivers a pre-set number of cells 
with pre-set characteristics. The pulse processor 
measures pulse width which is used to discriminates 
single cells from aggregates. 5 

Prior to sorting, sort accuracy was checked as 
follows. An aliquot of the pre-B cells was labelled with 
nuclear UV excitable dye Hoechst 33 342 (10^g/ml) 
for 10 minutes at 37° C for examination of cells under 
a fluorescence microscope. After washing and sus- 10 
pending the pre-B cells in PBS, a defined number of 
cells were sorted using a FACStarP' 113 based upon for- 
ward and side scatter gates to exclude dead cells 
which have lower forward and side scatter. The num- 
ber of cells sorted was verified under the fluores- 15 
cence microscope. Sorting efficiency was deter- 
mined to be greater than 97%. 



Example II 



20 



250 pre-B cells were sorted into 2.5 x 10 6 PBMC. 
The pre-B cells were stained with Hoechst 33 342 as 
in Example I and with monoclonal antibodies as de- 
scribed above. Sequential aliquots of cells were tak- 
en out of the staining solution just prior to analysis. 25 
The cells were washed twice, resuspended in PBS 
and were immediately analyzed on a FACStarP ,us cells 
sorter which had been cleaned as described above. 
Cells were sorted onto a microscope slide if they were 
FITC + (Le., FL1 + ), PE + (Le., FL2 + ) and PerCp" (Le., 30 
FL3"). Cells were sorted onto a new slide after every 
10 6 cells analyzed. Slides were examined under a UV 
fluorescence microscope. The presence of pre-B 
cells was verified on the slide by UV excitation. This 
experiment confirmed that the sort strategy worked. 35 

Example III 

0, 250, 2,500 and 25,000 pre-B cells were sorted 
into 2.5 x 10 8 PBMC. This represents a frequency of 40 
pre-B ceils of 10- 6 , 10 -5 and 10 -4 respectively. The 
pre-B cells were stained as in Example I with the first 
and second markers described above. The cells were 
washed twice, resuspended in PBS and were imme- 
diately analyzed on a FACStarP' us cells sorter which 45 
had been cleaned as described above. 

Referring to FIG.s 3H and I respectively, gates 
were established in a plot of FL2 versus FL1 and FL3 
versus FL1 for pre-B cells at a concentration of one 
in 1 CM PBMC. The gates, referred to as "R5" in H and so 
"R3" in I, were set so as to include the most likely 
areas were pre-B cells should occur. 

The R5 and R3 gates then were applied to the list 
mode data from the analysis, and a plot of SSC ver- 
sus FSC was constructed. Referring to FIG. 3A, the 55 
gate "R1 " was constructed in the scatter plot so as to 
include all of the pre-B cells. All three gates then were- 
applied to various aliquots. 



Referring to FIG. 3B and C respectively, PBMC 
were analyzed using anti-KLH as the marker without 
and with the R1, R3 and R5 gates. As can be seen, 
for ungated data numerous events are displayed in a 
plot of SSC versus FSC; however, when the gates are 
applied, only one event in 10 6 events analyzed is re- 
corded. Referring to FIG. 3D, when 0 pre-B cells are 
seeded into PBMC and 10 6 events are analyzed, only 
one event in 10 6 appears as background. 

Referring to FIG.s 3E, F and G respectively, 
gates R1 ; R3 and R5 were applied to list mode data re- 
corded from the analysis of 10-*, 10~ 5 and 10" 4 pre-B 
cells seeded into at least 10 6 PBMC. For 10" 6 seeded 
cells, 12-17 events were recovered in 4 replicate ex- 
periments. For 10" 5 seeded cells, 45-136 pre-B cells 
were recovered in 3 replicate experiments. For 10 -4 
seeded cells, 371-1291 pre-B cells were recovered in 
3 replicate experiments. As expected, the number of 
recovered cells increases by a factor of 10 as the 
number of seeded cells increases by a factor of 10. 

All publications and patents mentioned in this 
specification are indicative of the level of ordinary 
skill in the art to which this invention pertains. All pub- 
lications and patents are herein incorporated by ref- 
erence to the same extent as if each individual pub- 
lication or patent was specifically and individually in- 
dicated to be incorporated by reference. 



Claims 

1 . A method for the analysis of particles in a sample 
at frequencies of less than one in 10 4 comprising 
the steps of: 

(a) combining the particles with at least one 
first marker, wherein said first marker is spe- 
cific for the particle to be analyzed, and at 
least one second marker, wherein said sec- 
ond marker is specific for a majority of all 
other particles in the sample but not for the 
particle to be analyzed; 

(b) setting one or more fluorescence thresh- 
olds above the autofluorescence of the unlab- 
eled particles in the sample so as to exclude 
from analysis all particles that do not exceed 
that threshold; 

(c) setting a gate for each of the light scatter, 
first marker and second marker parameters 
measured so as to include the rare events to 
be detected and exclude all other events; and 

(d) analyzing the particles in the sample for 
on or more light scatter parameters and for 
the presence of the first maker and absence 
of the second marker by means of flow cyto- 
metry within the gates set by step (c). 

2. The method of claim 1 wherein the particles com- 
prise cells in a biological fluid selected from the 
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group consisting of blood, urine, lymph, bone 
marrow, ascites, cerebrospinal fluid and saliva. 

3. The method of claim 1 wherein the particles com- 
prise cells dis-aggregated from an organ or solid 5 
tissue. 

4. The method of claim 1 wherein the particles com- 
prise cells and the markers comprise immuno- 
fluorescence markers. 10 

5. The method of claim 4 wherein the number of first 
markers is two and each marker is labelled with 
a fluorescent dye having an emission wavelength 

that is distinguishable from the other. 15 

6. The method of claim 4 wherein the number of 
second markers is more than one and such 
marker is labelled with the same fluorescent dye, 

said dye having an emission wavelength that is 20 
distinguishable from the wavelength of the dye 
on the first marker. 

7. The method of claim 4 wherein the number of first V ..,if^» 
markers is two and the number of second mark- 25 r -s^zm 
ers is more than one wherein each of said mark- . * 
ers is labelled with a fluorescent dye, said dyes . 

on the first markers having an -emission wave- "s-.ts 
length that is distinguishable from the wave- 
length of the dye on the second marker(s) and 30 ' . . _ 
said dye on the second markers having an emis- ;.: 
sion wavelength that is distinguishable from the * 
wavelength of the dyes on the first markers. : . vt^s 

8. The method of claim 1 wherein step (c) is elimin- 35 ~ y- 
ated. 

9. The method of claim 1 wherein first marker is for 
a cancer cell and the second marker is CD34. 

40 

10. A method for the elimination of spurious-events 
from an n-dimensional datastream, wherein said 
data on each event is acquired by means of flow 
cytometry, comprising the steps of: 

(a) calculating the a priori probabilities of 45 
each event based upon the whole data file; 

(b) calculating the probability of each se- 
quence of 1 000 events based upon the a priori 
probabilities of step (a); 

(c) removing ten percent of the least probable so 
sequences as calculated in step (b); and 

(d) analyzing the resulting data for the pr s- 
ence of rare events. 
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